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An adiabatic approach is developed for the problem of boundary friction between two atomically 
smooth and incommensurate solid surfaces, separated by a monolayer of lubricant atoms. This 
method permits to consider very slow macroscopic motion of the parts in contact, separately from 
fast thermic motions of individual atoms. A characteristic "stick-slip" behavior of the tangential 
force on the contact is obtained within a simple ID model, relevant for the tip and sample system 
in friction force microscopy (FFM). This behavior reflects the specific mechanism of stress energy 
accumulation, through formation of long-living metastable states (defects) within the monoatomic 
lubricant layer, and their subsequent collapse with energy conversion into heat. This is similar to 
the dislocation mechanism of irreversible deformation in bulk solids. The peculiar feature predicted 
by the present theory is the twofold periodicity of "stick-slip" spikes with relative displacement: 
the shorter period ad (where a is the tip lattice periodicity and 5 the relative tip-sample lattice 
mismatch) relates to defect skips by one elementary cell, and the longer period a(l — 5) relates to 
defect annihilation or nucleation at the boundaries of contact area. 

PACS numbers: 07.79.Sp, 68.55.L, 81.40.Pq 



INTRODUCTION 



Atomic mechanisms underlying the phenomenon of 
friction attract now a still increasing interest of many 
investigators. Particularly, this is motivated by the need 
of better understanding the physical principles that gov- 
ern the image formation in friction force microscopy 
(FFM) [1]. Commonly, we call friction the process of 
energy dissipation during the relative displacement of 
two solid surfaces in contact [2]. The friction force 
F(x) = —dW(x)/dx, where W(x) is the work performed 
to reach the position x (Fig. 1), is always directed con- 
trary wise to the displacement. The empirical "Amon- 
tons' law" relating this force to the normal load N on 
contact: F — kN, where the friction coefficient k de- 
pends only on the nature of the contacting materials, is 
known for about three centuries [3]. The modern tribol- 
ogy much dealt with the linear relation between N and 
the effective contact area A, through deformations of mi- 
croscopic rugosities [4], so that F = S C A with a certain 
load independent constant S c (the Bowden- Tabor law). 

Nowadays, it is being progressively conceived that es- 
sential physics is already contained in the fundamental 
process of friction between two atomically smooth solid 
surfaces, not accompanied by any morphologic changes 
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FIG. 1: Friction force F vs applied load N at given velocity 
v of relative displacement of two bodies. 



(wearless friction). It is also well recognized that a rather 
necessary premise for such a regime is the presence of 
some intermediate (lubricant) layer of relatively inert 
atoms between the surfaces. Usually, this friction regime 
is referred to as boundary friction [2]. The most evident 
source of dissipation at wearless friction is the genera- 
tion, by microscopic interactions at the contact, of var- 
ious kinds of quasiparticles. Those should bear away a 
part of mechanical energy from the external drive and 
finally transform it into heat. 

To the date, a plenty of experimental material on this 
sort of friction is obtained. Particularly, a lateral force 
contrast was discovered in FMM scans along an atomi- 
cally smooth surface containing domains of different ma- 
terials [5, 6], thus showing the constant S c in the Bowden- 
Tabor law to be material dependent. Furthermore, in 
FFM experiments with the best resolution, the friction 
force vs time (or displacement) isn't constant but re- 
veals a sawtooth-like, or the so-called stick-slip, behavior 
[1]. This should indicate ocurrence of certain mechan- 
ical instabilities during the displacement, either within 
the contact itself or/and within the measurement appa- 
ratus (cantilevers, drives, etc.). The theory related to 
cantilever instabilities [7, 8] has its main disadvantage in 
that it doesn't account for the irreversible (microscopic) 
processes. The models for the contact instabilities, such 
as the independent oscillator model [9, 10] or the Frenkel- 
Kontorova model [11, 12], though being able to explain 
the energy losses, were limited so far to only direct solid- 
solid interactions (see also [13]), not mediated by the lu- 
bricant. For the processes of boundary friction, the most 
elaborated approach is that considering the sequence of 
"freezing-melting" transitions within lubricant layers of 
few atoms thickness as a source of stick-slip disconti- 
nuities [14]. However, the FFM conditions can already 
reach to ultrathin boundary layers, monoatomic or even 
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submonoatomic. In this situation, the interlayer viscos- 
ity and hence the melting transitions are most proba- 
bly excluded, nevertheless the indication of still present 
discontinuous friction forces (atomic contrast in FMM) 
poses a challenge for atomic friction theory. The treat- 
ment in this paper just concerns the contact instabilities 
between two solids separated by a monoatomic lubricant 
layer and define microscopic mechanisms for irreversible 
energy losses with atomic periodicity at slow relative dis- 
placement. To fulfill this program, it is of special impor- 
tance to choose an adequate calculation method. 

The most straightforward method, very popular now 
for modelling various processes on nanoscopic level, is 
that of molecular dynamics (MD), and it was also applied 
to boundary friction [15-17]. However, this method has 
one substantial limitation that was not yet addressed ex- 
plicitly. The time step for integrating the MD equations 
is typically of order of femtoseconds, then, in order to 
meet reasonable computing time requirements (no more 
than 10 4 - 10 6 steps), one has to consider the displace- 
ment velocities v not slower than ~ 10 3 cm/s to model 
the displacement by only few atomic periods. But in re- 
ality, even the faster process of atomic sliding of a free 
Lennard-Jones adsorbate layer on a single metal surface, 
subject to a uniform parallel force [18], is characterized 
by much slower velocities v ~ 0.1 cm/s. And the veloc- 
ities in STM, AFM, FFM, etc. techniques, scale from 
10~ 4 down to 10~ 7 cm/s, thus opening a 10 orders of 
magnitude abyss beyond the MD capacities. Meanwhile, 
just such slow adiabatic motions, as will be also seen 
further, arc characteristic for the observable friction phe- 
nomena. 

The remedy for this problem can be found in focusing 
the treatment on possible long living metastable states 
in the boundary layer, similar to those, well-known in 
physics of solids, as dislocations, domain walls, etc. Then 
the slow dynamics (SD) of such states can be considered 
separately from fast thermal motions of MD. Adiabati- 
cally, the SD part just follows the scenarios of contact 
instabilities of Refs. [9-13], while the MD part is ac- 
counted for in average, through the temperature depen- 
dent factor, triggering the transitions between different 
(meta)stable states. (This also can be meant as if the 
adiabatic atomic potentials were defined only within to 

Each transition is preceded by a certain period when 
the energy of a metastable state (m-state), and so the 
elastic force on the system, is growing with the macro- 
scopic displacement x, while the energy barrier h, sepa- 
rating the m-state from the nearest stable (s-) state, is 
decreasing (Fig. 2). The adiabaticity condition would 
require that the m-state lifetime, r m = T a exp(f3h), docs 
not change too rapidly (see a more precise criterion at 
the end of Sec. 3). When r TO decreases so much that it 
becomes comparable to the characteristic time t v = a/v 
of slow displacement by an atomic period a (though still 
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FIG. 2: Energy level scheme for transition between a 
metastable m-state to a stable s-state, through an excited 
s*-state. The excitation energy is e* and the barrier height 
h. 



at (3h ^> 1), the barrier can be overcome by a thermal 
fluctuation. Then, within a very short time ~ r a , the 
system passes from m to another state s*, highly excited 
by an energy e* 3> /3 _1 (see Fig. 2) over s. Next, within 
a certain relaxation time r r ~ 10 8 s, the energy differ- 
ence e* is released through emission of about (3e* 3> 1 
quasiparticles. 

Since r r is much shorter of the characteristic time 
(~ T m) of the inverse s* — >m transition, the latter pos- 
sibility can be safely excluded [19]. The considered m 
—*8* — > s process gives an elementary contribution to ir- 
reversible losses and can be compared to an individual 
spike in the stick-slip picture. The work by displacement 
Ax is written, accordingly to the 1st law of thermody- 
namics, as 

AW = AE + AQ (1) 

Between the transitions, the heat transfer AQ = and 
the force varies continuously (the stick stage) . Contrary- 
wise, on a transition, the entire energy change is released 
into heat: AE = —AQ , hence AW = 0, and a force 
discontinuity is generated (a slip). In this first attempt, 
the model is reduced to the simplified situation of ID in- 
terface between 2D "solids", however its extension to the 
realistic 2D (and possibly curved) interface can be also 
done on the same conceptual basis. 
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standard Lennard-Jones potential: 



/lj (y) = y~ 12 -y- 6 . 



(3) 




Tip ► 

• ••••••••• «o 

• • • •^•L • Lubricants 

oooooo ooooo oooo 
O O O O O C TO ooooo o o o o o 

•••••••••••••••• 

Sample 



FIG. 3: General view of the contact area. The dark circles 
represent tip atoms, the pale grey are sample atoms, and the 
empty circles are lubricant atoms. 



FORMULATION OF THE MODEL 

Let us consider a system (Fig. 3) consisting of two 2D 
"solid" arrays: conventionally called "tip" (t) and "sam- 
ple" (s), separated by an atomic chain of "lubricants" 
(0 at some distance d. We assume the same triangular 
lattice for both solids, but admit a small mismatch 5 be- 
tween the lattice parameters. The i-side of the contact 
consists of L elementary cells (L + 1 atoms), while the 
s-side is unlimited, and the Z-chain in contact with t con- 
sists of L atoms. Without loss of generality, we can put 
the i-lattice parameter equal a and the s-lattice parame- 
ter <z(l - 5), S < 1. 

These lattice structures are supposed not to vary at 
bringing t and s into (Z-mediated) contact and during 
their posterior mutual displacement by x. Thus we as- 
sume the solid-solid links infinitely rigid compared to the 
solid-lubricant and lubricant-lubricant links. Hence, the 
only variable part in the full system energy is that related 
to the Z-subsystem: 



n=l 



U (r n ,x) 



(2) 



where U(r n ,x) is the "mean-field" potential on the n- 
th lubricant with 2D radius-vector r„ for given t- s dis- 
placement x, and V (r n . n :) is the l-l interaction (r nj „/ = 
r n — >v). Strictly speaking, the energy, Eq. (2), is also 
a function of all r„, but if we only admit them to be- 
long to stable (or metastable) Z-configurations at a fixed 
distance d, the displacement x remains the single rele- 
vant parameter. The specifics of the above formulation 
consists, firstly, in taking into account two atomic period- 
icities (generally incommensurate) at once and, secondly, 
in the bilateral restriction of lubricant layer between two 
surfaces. Qualitative difference of this situation from the 
known microscopic models of a lubricant layer over sin- 
gle solid surface with single periodicity [9-13, 17] was 
already recognized in literature [20]. For simplicity, we 
will construct all the interactions in Eq. (1) using the 
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U (r n ,x) = s 
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where Eq is the adhesion energy, r m and r m i are re- 
spectively the coordinates of t- and s-atoms, r t and r s 
the corresponding t-l and s-l equilibrium distances, and 
x the vector of displacement of the t-. Thus, both 
the t-l and s-l coupling energies equal Eq, serving as 
the energy scale. At least, the l-l interaction is taken 
V (r n , n <) = ,ge /Lj (r n ^/r t ), where g is the coupling con- 
stant and n the l-l equilibrium distance. 



MEAN-FIELD POTENTIAL AND METASTABLE 
STATES 

One can reasonably consider the l-l interaction the 
weakest one in the above system and, at the first step, 
put g — (this is a good approximation unless the lu- 
bricants approach each other too closely compared with 
n). Then the equilibrium states will correspond to vari- 
ous distributions of L lubricants over the minima of the 
mean-field potential produced by the t- and s-lattices. In 
what follows we suppose that the infinite sample moves 
with velocity v in a;-direction with respect to the fixed 
tip, and the tip-sample distance d is constant. We de- 
fine the n-th elementary cell of the boundary layer as the 
rectangle limited by the n-th and (n + l)-th tip atoms 
in x-direction, and d wide in y-direction, and the entire 
contact area corresponds to 1 < n < L (inset to Fig. 4). 

Let the initial tip-sample relative position x = at 
the instant t = be such that the s-atoms are located 
symmetrically in the zeroth cell. In what follows, we 
consider the potential energy U n of a lubricant in n-th cell 
as a function of x n only (as shown in Fig. 4), supposing 
y n be always adjusted to the relative minimum of this 
energy at given x n (i.e., supposing the lubricants only 
to move along energy "valleys" in 2D cells). Because of 
the small mismatch <5, the static potential relief U n {x n ) 
(0 < x n < a) slowly varies (in general, incommensurably) 
from cell to cell along the interface so that the cells aren't 
equivalent. If we neglect the boundary effects (as is the 
case in what follows), the potential relief within n-th cell 
is fully determined by the "phase parameter" p n — x — 
n<5(mod(l — 5)). 

As seen from Fig. 4, besides the general raising or low- 
ering of minima with p (elastic strains), there also exist 
some critical values: p c \{~ 0.22 in that particular case) 
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FIG. 4: Mean-field potential profiles for a lubricant atom in 
the unit cell of a contact, at different values of the phase 
variable p (marks at each curve). The mismatch parameter 
S = 0.05, the interfacial separation d = 1.12a. Inset: the 
phase variable p is defined as the distance along x-coordinate 
between the midpoints of nearest neoghbor t- and Z-atoms 
within the cell. 
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and p C 2 = 1 — S — pd , when the initial single minimum 
£A°) splits into two. In the splitted potential, the lower 
minimum is separated from the upper one by a 
maximum U^ 3 ' (a saddle point in the 2D cell). The slow 
relative t-s displacement will result in that the potential 
relief as if "moves" adiabatically along x with the "phase 
velocity" v/S. 

If all the lubricants always rest in the lower minima 
the distances between them never become too short 
in this process and we can still neglect the l-l interactions. 
Then no irreversible losses and, hence, no friction occur in 
the system (the force on external drive F — ~dE(x)/dx, 
though non-zero, results fully elastic). However, such an 
/-configuration, even at moderate L (note, however, that 
a contact area in FMM typically includes ~ 10 3 or more el- 
ementary cells), cannot persist for any macroscopic time 
by the thermodynamic reasons. 

Indeed, with growing displacement, the cell poten- 
tial splits and then the barriers h\ — — and 
/i2 = — begin to grow and eventually get prac- 
tically impenetrable. Then there is a finite probability 
w > that at this stage a lubricant atom leaves "cap- 
tured" in the upper, metastable minimum U<- 2 \ For in- 
stance, in the model of linearly growing barriers with 
time (see Appendix A), this probability is simply given 
by the Fermi distribution function for the two-level sys- 
tem: 

w(AE) = (c^ + iy 1 (4) 



FIG. 5: Relative positions of two interacting lubricant atoms, 
n and n + 1, against the mean-field potential (upper part) 
and the valley of steepest descent in xy-plane (lower part). 
The mean-field potential is that of Fig. 4 and the interaction 
constant g = 1/120. The phase p n is close to the transition 
value pt . 



being AE the energy difference between the two minima. 
If all L lubricants are initially in states, the time 
necessary for at least one of them pass to is about 
l/(Lwv), which corresponds to the system displacement 
by only l/(Lw) atomic periods while Lw can easily ex- 
ceed unity. 

(2) 

Now let the n-th lubricant leave "captured" at Uk , 
the upper minimum of n-th cell, while its neighbor in 
(n + l)-th cell stays at U^h, and follow the system 
evolution with displacement (Fig. 5). As the distance 
Tn,n+i gets shorter, this particular l-l interaction cannot 
be more neglected. Such a pair of strongly interacting 
lubricants just realizes a metastable state for the consid- 
ered system. The Lennard- Jones repulsion will slightly 
shift the 7i, (n + 1) pair from the mean-field minima Un 

and towards the saddle points Un^ and U^+i- If 

the interaction constant g isn't too small (for the par- 
ticular case considered in Sec. 4, g = 1/120 is already 
sufficient), such repulsion makes the energy barrier be- 
tween this m-state and the s-state to turn from growing 
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to lowering. The displaced equilibrium positions U n 
and U n \\ are denned by the compensation between the 
I- 1 interaction and the crystalline field forces along the 
valleys: 

A n U n ■ (A n U n + A n U rhn+1 ) = 0, (5) 
A n +iU n +i ■ (A n+ iU n +i + A n +iU n ,n+i) = 0, 

(A n stands for the 2D gradient along r n ). Notably, the 
numeric results (see Fig. 5 and Sec. 4) show that the 



energy difference U^ + \ — U^li for the (n + l)-th lubri- 

cant increases much faster than Un 
one. Hence, as the repulsion grows with x, the barrier 
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(2') 



U n 2 ^ for n-th 



l n+l 



Unli ~ U n x + { decreases, while h£' = U n °> - U n ' 
remains increasing, and the (n + l)-th atom will be fi- 
nally pushed out by a thermal fluctuation from the uj^\ 

(2) 

position to Then the distance r n) „+i will suddenly 

r< 

increase and therefore the n-th atom will relax from U n ' 

(2) 

to U n ■ The overall energy gain 

Ae - U^ + U^ - IJW - U {2) 



r(i') 



(2) 



r(3) 



r(2') 



(2') 



n+1 



(not to be confused with the energy difference AE in 
Eq. (4)) is released through emission of quasiparticles. 
Supposing isotermic conditions, this energy is eventually 
transferred from the external drive to the thermostate. 

Further on, the n-th lubricant stays always at U n (un- 
til it returns to the unsplitted position Un^ at p n — * p C 2) 
while the metastable pair is now formed by the (n + 1)- 
th and (n + 2)-th lubricants (that is, skipped by one 
cell) and so on. For the hierarchy of times specific for 
our system, the transition time t tr (related to the corre- 
sponding phase value pt r through p tr = vt tr ) has a very 
small statistical dispersion around its mean value (tt r ) 
(see Appendix B). Also the mismatch parameter 6 can 
be chosen so that either the "capture" and "skipping" 
processes affect only single metastable pair at a moment, 
not its nearest neighbors, thus avoiding possible compli- 
cations due to interacting fluctuations. 

Once randomly initiated, each skipping process de- 
scribed above provides a rather regular and stable gener- 
ator of energy losses within a "phase domain" of length 
Id ~ a(l — 5)/ 5 in the Z-chain. Along such a domain, the 
phase p changes from to 1 — (5, and the lubricants are in 
the U^ position at all the sites withp„ < (ptr) — v (Ur), 
and in U^ position at p n > {ptr) (for more details see in 
6). The spikes in the friction force from subsequent skips 
have approximate periodicity rj Sa/v in time, or w Sa in 
displacement. Since this is much smaller of atomic peri- 
odicity, it may be difficult to resolve such fine structure 
in actual experiments. 

But the atomic periodicity itself can result from an- 
other source. The described generators move along the 
chain with phase velocity v p = v/S and, at sufficient dis- 
tances from the beginning of chain (n = 1), they exist 




FIG. 6: Domain structure of the lubricant chain formed at 
displacement of the contacting solids. The "spin" representa- 
tion is chosen, where the spin-up corresponds to an atom in 
upper well, spin-down to that in lower well, and spinless is an 
atom in unsplitted well. 



within each domain, Id apart. The friction force, aver- 
aged over the (5-periods, is (F) s ks n g As/(5a), where n g 
is the total number of generators in the chain. When a 
generator reaches the chain end (n — L), it disappears 
(like a dislocation when it reaches the sample surface) 
and n g decreases by unity. Hence a jump to down occurs 
in the average force (F) s . Otherwise, when a new gen- 
erator appears near the chain beginning, (F) s increases. 
These events have the periodicity w lg/Vp — 0,(1 — 5)/v 
in time, that is just the sample cell periodicity a(l — S) 
in displacement. 

The averaged friction force (F) over many atomic pe- 
riods results simply (F) ~ LAe, that is proportional to 
the "contact area" L with the Bowden- Tabor coefficient 
Ae, which depends both on the material parameters and 
the interfacial distance d (related to the pressure). 

To conclude this section, we obtain a quantitative cri- 
terion for adiabaticity in our system. The limitation on 
the displacement velocity v can be extracted from the 
evident condition that the lifetime r m of the metastable 
state (especially when it becomes comparable to t v ) 
doesn't changes appreciably within the microscopic time 
~ il^ 1 (atomic oscillation period): dT m /dt <C Or m , or 
[3dh/dt <C Since dh/dt = vdh/dx and the barrier 
height changes from its maximum value h ma x to zero 
within a fraction q <C 1 of atomic period, the adiabatic- 
ity criterion is given by: 



V <C V cr 



(6) 



The numerator in the last equation is of the order of 
sound velocity while the denominator (3h max /q can be as 
high as ~ 10 2 . Thus, the adiabaticity must be definitely 
violated (with all possible complications, as appearance 
of shock waves, thermal instabilities, etc.) at v above 
some m/s, as is the case for the most of MD simulations. 
However, it is well assured at actual FFM velocities, men- 
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Displacement (cell units) 



FIG. 7: Microscopic tangential force on the contact vs sample 
displacement. Each fine structure spike corresponds to a skip 
of the boundary between "spin-up" and "spin-down" areas of 
Fig. 6 by one cell period, while the "envelope sawtooth" with 
the periodicity of sample lattice (0.95a) corresponds to the 
disappearings of domain boundaries at reaching the sample 
end. 

tioned above. 

NUMERICAL RESULTS 

We choose the equilibrium distances for Lennard-Jones 
interactions r t — 0.8a, r s = 0.78a, r; = 0.9a, and the 
mismatch parameter 5 — 0.05. The energy parameter £o 
is chosen 10 meV and the cell parameter a — 0.3 nm. 
At least, the interfacial distance is chosen d = 1.12a. 
This value is slightly below the equilibrium separation 
d e = 1.24a for a cell at the displacement phase p = 
and corresponds to the normal force on tip N — 17.2 nN. 
The calculated profiles of mean-field potential for differ- 
ent values of the phase parameter p are shown in Fig. 4. 
Next, the equilibrium equations, Eq. (5), are resolved to 
give the values of barriers and hP^ as functions of p. 
The total length of the system was taken L = 20 and its 
initial state corresponds to one lubricant at J/W (or U^) 
state in each cell. The displacement velocity v = 30 nm/s 
corresponds to the characteristic time of displacement by 
cell period t v — 10~ 2 s, and the elementary time step is 
At = 10~ 5 s, that is 10 orders of magnitude slower of 
the respective MD times. The evolution of total energy 
E(x) is modelled, using the mean-field and interaction 
energies (at displaced positions for metastable pairs) at 
phase values p n which vary with x = vt, and the account 
for capturing and skipping events with probabilities is de- 
scribed in Appendices A,B. The temperature is set T = 1 
meV w 11 K. 

The resulting plot of the friction force F{x) is pre- 



sented in Fig. 7. After some initial period, when the force 
is almost elastic, the stick-slip behavior with two types 
of periodicity establishes and the mean level of friction 
force corresponds to a rather high "friction coefficient" 
k « 0.6. A more detailed study of the present model, 
for different values of pressure, temperature and velocity 
parameters, will be given elsewhere. It is also of interest, 
that using the above parameter values and f2 = 10 12 s -1 
in the adiabaticity criterion, Eq. (6), we get the criti- 
cal velocity v cr ~ 0.3 m/s, in good accordance with our 
previous reasoning. 

CONCLUDING REMARKS 

The above analysis gives only a simple illustration for 
possible irreversible processes within the adiabatic pic- 
ture of boundary friction. Other types of the system 
behavior can be obtained at varying the material and 
external parameters. In particular, with growing separa- 
tion d between surfaces (lower or even negative pressure) 
the potential profiles get lower and the stick-slip behavior 
is expected to disappear. Contrarywise, at closer sepa- 
ration, the stage of inverse barrier evolution ceases to 
appear for lubricant pairs, and irreversible transitions 
(much sharper) can occur already in close clusters of 
three or more lubricants, however this would greatly com- 
plicate the treatment. A special regime can be obtained 
for the case of two fully commensurate solids ((5 = 0). In 
this case, there is no preference between the two splittcd 
mean- field states in a cell, and the "phase domain" be- 
comes infinitely long. All the skips would occur simulta- 
neously over the whole contact area, at the moment when 
the transition phase value ptr is reached at all the cells. 
This will generate very high and short spikes of friction 
force with single periodicity of the lattice. Note in con- 
clusion that passing from the above considered ID to the 
realistic 2D situation will complify the model, not only 
by the extension of arrays for computation, but also due 
to introduction of several types of different metastable 
states. 

APPENDIX A 

To model the initial evolution stage of potential in our 
system, let us consider a particle in the double-well sys- 
tem (Fig. 8) where the potential barriers for each well 
hi 2 vary adiabatically with time at t > 0. 

Then the probability w(t) to find the particle in the 
upper well at the moment t obey the linear differential 
equation with the coefficients slowly depending on time 
(through hi, 2 (<)): 

fT 1 ^ = - (e~^ + e"^ 2 ) w + e-^. (Al) 
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FIG. 8: Double-well potential for an isolated lubricant atom. 



The solution to Eq. (Al), satisfying the initial condition 
w(0) = (that is, at t — the particle is with certainty 
in the lower well), is given by the evident expressions: 



where 



-0hl(t') +e -/3hi(?) 



dt'. 



(7) 



(A2) 



Now let us adopt, for simplicity, the model linear law for 
the growth of the barriers: 



hi (t) = AE + h 2 (t) , h 2 (t) = ut. 



(A3) 



As seen from Fig. 4, Eq. (A. 3) is a rather plausible ap- 
proximation. Then the explicit integration in Eq. (A. 2) 
gives the saught probability: 



w(t) = 



1 — exp [— f2 (l — e~ 



0Ut\ 



IH 



-3AE , 1 



(A4) 



Since the rate of barrier growth is u = vh max /(aq), we 
conclude from comparison with Eq. (6) that the dimen- 
sionless parameter Cl//3u is just the adiabaticity ratio 
v cr jv and, for typical FFM conditions, its value is enor- 
mous: ~ 10 8 . Thus the exponential in the numerator of 
Eq. (A. 4) rapidly vanishes and we arrive at the result of 
Eq. (4). 

APPENDIX B 

Let us turn to the stage of reversed evolution (see the 
paragraph before Eq. (5)), and let the potential barrier 
(the lower one in Fig. 8) adiabatically decrease in time as 
h(t) = —ut, — oo < t < 0. If at the moment t the particle 



FIG. 9: Probability distribution p, normalized to its maxi- 
mum value pmax , as a function of dimensionless variable (5ut 
(see text) at different values of the parameter v c /v (indicated 
at each maximum). The peak is quite narrow, its position 
varies very slightly over two orders of magnitude variation of 
v c /v, while its form is just invariable. 



is in the well, the probability that it escapes within an 
infinitesimal time interval dt is dw(t) = fle^ ut dt. Then 
the probability distribution p(t) for the transition to oc- 
cur at the moment t is the product of dw{i)/dt times the 
probability to survive within the well during the period 
from — oo to t: 



P(t) 



dw 
Ut 



exp 



dw {t') 



exp 



']3u 



But 



j3ut 



The function p(i), shown in Fig. 9 for different val- 
ues of u, sharply peaks at t* = ln(/3u/fi)/ j3u. Since the 
parameter VL/ (3u — v cr /v is huge [21], then: i) the distri- 
bution p(t) is very narrow and ii) its form and the critical 
barrier value h* = h(t*) — /3 _1 ln(v cr /v) depend on the 
displacement velocity v only very weakly. 
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